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I . INTRODUCTION 

The past five years have seen a revival of theoretical interest in 
predissociation phenomena, coupled with increasing sophistication in 
experimental technique [1] and a continuing search for related resonance 
effects in the molecular scattering field [2-5]. Given the established 
characterisations of predissociation types [6-7] and selection rules 
[8-9], and the accepted value of predissociation measurements in 
locating the dissociation limit [6,10], the object of recent theoretical 
work has been to underline further information available from individual 
line width (or lifetime) measurements. This information relates 
primarily to the forms of the relevant potential energy curves. 

The purpose of this review is to summarise our present under- 
standing of predissociation by rotation (type III [6]) and by curve 
crossing (type I [6]) in diatomic molecules. The material therefore falls 
into two parts ) each of which includes a formal quantum mechanical intro- 
duction, which may also serve to emphasise the connection between 
spectroscopic measurements and molecular scattering data. Subsequent 
sections cover the available computational and analytical techniques, 
followed by a detailed summary of conclusions. An attempt has been 
made to include all line width or lifetime calculations for specific 
systems published before going to press (May 1973). 

The first general conclusion is that the pattern of predissociation 
line widths provides a highly sensitive yardstick for the determination 
of unknown potential curves. Secondly, the computation of such a 
pattern for given potential curves is now a matter of routine, unless 
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the predissociation happens to occur from an upper adiabatic (avoided 
curve- crossing) potential curve. As an alternative to exact 
numerical computation one may expect to apply available relatively 
rapid analytical (semi-classical) estimates within an accuracy of 10%. 
These analytical formulae also to provide physical insight into the 
details of the predissociation pattern, to the extent that a direct 
inversion procedure has been developed for determination of the 
repulsive potential curves responsible for type I predissociations. 


II. PREDISSOCIATION BY ROTATION 

The relation between predissociation by rotation [6] and the shape 
or orbiting resonances of collision theory [2-5, 11] is by now well 
established. Both are due'to the existence of a quasi-bound state 
separated from the continuum by a potential barrier, as depicted in 
figure 2.1, but the two types of experiment are, in fact, complementary 
since the scattering resonances most readily detected experimentally 
correspond to lines which would appear indiscemibly diffuse in the 
spectrum. 

Much of the relevant theory is to be found either in the 
scattering literature, or at least couched in the language of scattering 
theory. Advances have been made by both computational [14-21] and 
semi- classical methods [13, 12-32]. Some computational workers have 
used the resonant jump in the scattering phase shift [14-17], or the 
closely associated peak in the collision delay time (see below) [20, 
33-35] to determine the positions and energy widths of the quasi-bound 
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states. Others have pointed to the spectroscopically important peak 
in the relative amplitude of the wavefunction inside and outside the 
barrier region [14, 18-20]; however the intimate connection [12] 
between this characteristic and the behaviour of the phase shift has 
received relatively little emphasis in the molecular literature. Again 
in the semi-classical field attention has concentrated on the behaviour 
of the phase shift [24-27], or equivalently (see below) on the location 
of complex energy poles in the scattering matrix [28-30], but the 
behaviour of the amplitude ratio has been largely ignored. 

Our first purpose will therefore be to demonstrate the necessary 
connection between resonant characteristics of the phase shift, q(E) > 
the collision delay time "^(E) and the amplitude ration, A(E) . This 
general discussion will involve the introduction of complex energy 
levels 


E = E - ir /2 , T > 0 (2.1) 

n n n 

and lead to Breit-Wigner parameterisations of 0(E) , t,(E) , and 

d 

A(E) . We shall then specialise to a semi-classical (JWKB) descrip- 
tion in order to relate the reasonance positions, E^ , and widths, 

T n , to the form of the potential function. Finally, the practical 
applications of the theory will be discussed. 

^2.1 Resonance characteristics 

Given that the quasi-bound state with angular momentum J , 
supported by a potential, V(r) , lies in the continuum it is convenient 
to take the energy zero at the dissociation limit, Vi 00 ) = 0 , and to 
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define 

k 2 = 2mE/fi 2 
U(r) = 2mV(r)/^ 2 . 


( 2 . 2 ) 


The wavefunction must then satisfy 


~ + k 2 - U(r) 
dr 


J(J + 1) 

2 

r 


X EJ ( r ) ~ 0 , 


(2.3) 


subject to the boundary condition x E j(0) = 0 . The most convenient 
normalisation is to a delta function of energy [36] ; 

C X EJ (r) X E'J (r)dr = ^ " E ') * ( 2 * 4 ) 

This fixes the asymptotic amplitude of the wavefunction as, 

1/2 

sin[kr - Jtt/ 2 + lj^(E)] , (2.5) 

where 1j(E) is the phase shift, the term -Jtt/ 2 being included in 
(2.5) so that rij(E) =0 in the absence of any distortion potential 
(V(r) = 0) . Closely associated with tlj(E) is the collision delay 
time function [33-35] 

3n T 

x d (E/)= 2ft ~~ , (2.6) 

the resonance behaviour of which is more conveniently characterised 
than that of the phase shift [20]. 

The significance of the energy normalisation given by (2.4), is 
that it ensures unit density of states [37]. Hence the differential 


X E,J (r) 


2m 

Tffi 2 k 
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oscillator strength for aspectroscopic transition from a bound level 
^(r) of some other electronic states varies as 

v|/ x b (r) M(r) x E>J (r)dr| Z , (2.7) 

where M(r) denotes the electronic transition moment. Thus, with the 
external amplitude fixed by (2.5), a peak in the internal to external 
amplitude ratio leads to strongly enhanced overlap with (r) > and 
hence to a peak in the spectrum. 

The necessary connection between the behaviour of this amplitude 
ratio and that of the phase shift is however most readily demonstrated 
by the introduction of a generalised solution of (2.3), 

renormalised at the origin such that [12] 


lim 
r -* 0 


-J-l 


X kJ^ r) 


1 , 


( 2 . 8 ) 


and capable of extension to complex values of k , 


k = k' + ik" , 


(2.9) 


and hence according to (2.2) to the complex energy levels given by 
(2.1) , with 

E = (k' 2 - k" 2 )fi 2 / 2m 
n n n 


T = -2ik’ k"ii 2 /m . 
n n n 


( 2 . 10 ) 


It may be shown [12] that this generalised solution may be expressed in 

the asymptotic region, r -*■ 00 , in terms of two Jost functions f (±k) , 

J 
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r -> oo 

x k j (r) ~ 2 l [ f j< k ) eikr -fj(-k)e~ ikr ] , ( 2 . 11 ) 

where for complex k 

fj(-k) = [fj(k*)]* , (2.12) 

and fj(k) is analytic in k for Im k < 0 . It follows that for 
real k (and hence real E ) . 

r *■>* oo 

X kJ ( r ) ~ | f j(k) | sin(kr - .Jtt/ 2 + rjj (E) ] , (2.13) 

where 


n J (E) = arg fj(k) + Jtt/2 , (2.14) 

in agreement, apart from a change in normalisation with equation (2.5). 
This shows that both the amplitude ratio (which varies by virtue of the 
new normalisation condition (2.8), as |fj(k)| "S and the phase shift 
depend on the behaviour of the function f,(k) . 

J 

Particular attention attaches to complex roots of the f (±k) , 

vJ 

of which two types may be recognised. The simplest are zeros of f,(k) 

*J 

on the negative imaginary axis, k = - iY vJ say, so that according to 

( 2 . 8 ) 


*u (r) 


A e 


-V 


(2.15) 


These therefore correspond to the bound states with real negative 

2 2 

energies E vJ = -Y v j & /2m . The quasi-bound states on the other hand 
are associated with roots of f j(~k) , given by (2.9) with k' > 0 , 
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k" < 0 . Hence r given by (2.10) is positive. For roots of this 

type, ¥ (r) has purely outgoing characteristics 

K.J 

r * “ 1 ikr 

\,< r > ' k f J< k > e • <2 - 16) 

and as might be expected for a quasi-bound state, the imaginary 
component in the energy given by (2.1) ensures an exponential decrease 
in the time evolution factor 

|exp(-iEt/iO 1^ = exp (-T^t/ii) , (2.17) 

with time constant 

T = &/T . (2.18) 

n n 

Complex energy (and hence nonphysical) solutions of this type are 
termed Siegert [38] states, while the energies themselves are termed 
poles in the S matrix, because with the elements of S defined by 
the ratios of the outgoing to the incoming amplitudes of the wavef unction, 

Sj = f j(k) /f a (-k) , (2.19) 

in the present single channel problem. Hence, has a pole at every 
zero of fj(-k) . 

A linear expansion in the neighbourhood of such a pole, 

f T (k) = [f_(-k)]* = F(E)(E - E - ir /2) , (2.20) 

J J n n 

where F(E) may vary slowly with E may now be used to demonstrate 
the characteristic resonance behaviour in the physical (real energy) 
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solutions of (2.3). Clearly according to (2.14) and (2.20), with 
A(E,J) used to denote the amplitude ratio 

|A(E,J)| 2 - |F(E)| _ ^(E - E n ) 2 + r 2 / 4] _1 (2. 21) 

H,(E) = n°(E) + arctan[T /2(E - E) ] , (2.22) 

J J n n 

where 

r)j(E) = argF(E) + Jtt/2 , (2.23) 


while the collision delay time defined by (2.6) becomes 

30j 

x d (E,J) = -jg- + 


r -ft 

n 


(E - E ) 2 + T 2 /4 
n n 


(2.24) 


These are the standard Breit-Wigner parametrisations of A(E,J) , 

Tlj(E) and x d (E,J) . 

Clearly, if F^ is sufficiently small that variations in F(E) 
may be ignored over the resonance region, both the spectroscopic line 
shape and the collision delay time functions follow the same Lorentzian 
form. The line width parameter F^ may therefore be taken as the full 
width at half height of either function. may also be deduced from 
the maximum resonance increment, 455/r^ > in the time delay function, or 
in a purely analytical treatment (see below) from the imaginary part 
of the pole in the S matrix. 

Possible distortions in the line shape and changes in the peak 
position as the line width increases have recently been noted [20]. 
Since such distortions to A(E,J) and T d (E,J) depend on | F ( E) | 
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and ang F(E) respectively, they may be expected to differ in kind. 

Hence in any comparison with experiment A(E,J) must be employed in 

spectroscopic, and T^(E,J) (or rjj CE) ) in scattering applications. 

These general considerations therefore underline the intimate 

connection between the resonance behaviour of Ht(E) , x ,(E.J) and 

J d 

A(E,J) , and indicate the origin of the Breit-Wigner parameterisations 
given by (2.21)-(2.24) in the sharp resonance limit. The next step is 


to relate the resonance positions E and widths T to the form of the 

n n 


potential function V(r) . 


§ 2.2 Semi-classical theory 

The semi-classical (JWKB) method of solution of equation (2.3) 
relies on establishing connections between solutions of the form [39, 42] 


1/9 f i/^ k T (r)dr -i/ r k (r)dr 

X EJ (r) = [kj( r )r 1/Z X' e + X"e X 


(2.23) 


where [40] 

kj(r) = k 2 - U(r) - (J + \ ) 2 /r 2 , (2.26) 

which are valid in the regions a « r << b , r >> c in figure 2.2, 
where 



kj(r) 


« 1 . 


(2.27) 


Such connections are based on more accurate model solutions of (2.3) 
in the non-classical regions around the turning points a, b, and c 
where k^(r) = 0 . 
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Changes in will be followed by setting (X', X") in (2.25) equal 

to (A', A"), (B ’ , B") or (C', C") according to the choice of phase 

reference point x as a, b, or c ; single and double primes are used 

to denote outgoing and incoming motion respectively, as illustrated in 
figure 2.2. This means since the same function may be referred either 

to a or to b in the region a « r « b that 


ia. 


' B' 


B" _ 



-xo, 


r a- i 

A” 


(2.28) 


where 


a J ‘ Ja k J (r)dr 


(2.29) 


The significance of the superscript zero will become apparent later. 

The correct behaviour of the wavefunction at an isolated left- 

hand turning point, a , now requires that \p„ T (r) should behave as 

.hJ 

the asymptotic form of an Airy function [41], which is the exact 
solution of (2.3) for a linear approximation to the potential, decreasing 
exponentially into the non-classical region r < a . The appropriate 
form is [39, 42] 


X EJ (r) = A[k J (r)] 1/2 sin(/^ kj(r)dr + tt/4) . (2.30) 

Hence by comparison with (2.25) (with (A’, A") in place of (X', X")), 


A' 


j Ae - 11,/4 


A" - )As i,T/4 


(2.31) 
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Equally, at energies below the dissociation limit for which tjj p j(r) 
must decrease exponentially for r > b , it may be shown [42] that 

XgjCr) = BCkjCr)] -172 sin(/^ k a (r)dr + tt/4) , (2.32) 

so that 


B' = 


1 » iTr / 4 
2 Be 


B" = 


1 -n -iTT/4 

2 Be 


(2.33) 


Consistency between (2.28), (2.31) and (2.33) requires that the energy 
shall be chosen to satisfy the Bohr quantisation condition 


a ° = £ k j (r)dr = (n + i ) ' 


(2.34) 


for a bound state. A semi-classical normalisation of the wavefunction 
[42] then yields values for A and B in (2.30) and (2.32), 


n 


A = (-1) B = 


2mhoJ 

vfl 2 


1/2 


(2.35) 


where /fioj is the local energy level spacing, given by the quadrature 

■ I! = [ Ja 

Analysis of the quasi-bound states requires a proper treatment of 
the barrier region, in place of a simple boundary condition at b . 
Different models for the barrier yield different relations between 
(B', B") and (C ’ , C") in figure 2.2 [22-30], but any connection 
formula is subject to the equations of flux conservation 

| C ' | 2 - | C" I 2 = | B T | 2 - | B” | 2 , 


(2.37) 
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and of symmetry under time reversal, so that the substitutions 

(C 1 , C") -> (C"*, C'*) imply (B', B") -*• (B"*, B ' *) . 

(2.38) 

Hence we may express the general connection in the form 




(2.39) 


and give specific expressions for the parameters y, 0 and <j) for 
different barrier models in equations (2 . 59 )- (2 . 76) below. Solution 
of (2.39) when B' = 1, C" = 0 gives the following expressions for the 
barrier transmission and reflection probabilities 



tanh^y 


(2.40) 


It now remains merely to combine equations (2.28), (2.31) and (2.39) 
in order both to determine the phase shift and to express the external 
amplitudes (C 1 , C") in terms of the internal coefficient A ; 

coshy e * , -i sinhy e 

• . ~ “16 , ~ i<j> 

i sxnhy e , coshy e 
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[e^cos ctj + i e r sin a^] e v 


■ -Y • 

~t 1 r* n i 


£± /#e - i - tt/2) 


Y 

[e cos 


. -y . , -6./2)ce - $ - tt/2) 

a T - i e 'sin a T ] e ^ A r 


— A (cosh2Y + sinh2Y cos 2a^) 


it$ 


-15 


(2.41) 


where 


- 1 a + e ) 


6 = arctan(e tan a) + ~ (0 - $) - tt/4 


(2.42) 


This means that, with |c'| chosen to satisfy the asymptotic 
normalisation given by (2.5) 

1/2 


^EJ (r) 


2m 

•nh 2 k 


sin[kr - Jtt/ 2 + n°(E) + nj res) (E)] , 


(2.43) 


where 


n>) 


(res) 


lim {/* k(r)dr - kr + (J + )tt/2 + j (9 - $)} (2.44) 

r oo “ ^ 


"2y 


Pj ' (E) = arctan(e 1 tan a) , 


(2.45) 


while in the internal region, a « r « b , 


XT-,,(r) = [cosh2y + sinh2y cos 2a] 


- 1/2 


2m 


7* k(r) 


1/2 


sin (/ k(r)dr + 

cl 


tt/4) 


(2.46) 
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This differs from the bound state solution given by (2.30) and (2.35) 
by the factor 

A(E,J) = [fuo(cosh2y + sinh2y cos2otj)] _ (2.47) 

2 

Hence |A(E,J) | governs the spectroscopic line shape. 

Clearly, at least for large y when according to (2.40) the barrier 
transmission probability is low, resonances in q (E) and A T (E) occur 
at energies E j such that 

aj = f h a k J (r)dr - | (0 + <j>) = (n + | )tt . (2.48) 

This may be regarded as the quantisation condition (2.34) for a bound 
level subject to a small (see below) level shift due to the phase terms 
(0 + ^) ; as a rough estimate based on a linear expansion for a° and 
neglect of any energy variation in (0 + (j)) , 

AE nJ * E „J - E L “ <~ 6 „J + • < 2 - 49 > 

where E°j is the unperturbed level given by (2.34) and /fito is defined 
by (2.36). The subscripts imply that 0 and are evaluated at the 
resonance point. E^j may also be obtained directly (and more accurately) 
from (2.48). 

It may be seen with the help of (2.45) and (2.47), to the extent 

that energy variations in y are unimportant over the resonance region, 

first that the resonance part of collision delay time T,(E) and the 

d 

line shape |a(E,J)| have the same functionality 
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res - „ „ S.n (res ^ 

T d (E,J) - 2B 


= 27ft! 1 A(E, J) 


= (2tt/uj) [ cosh2y _ + sinh2y cos 2a ] ^ , 

nJ nJ J 

(2.50) 


and secondly that the total integrated intensity from a single line is 


equal to unity 


E T +ihw 
nJ 2 


E '-kw 
nJ 2 


A(E,J)| 2 dE * ^ 


3E 


3a, 


''(n+l)7T 

nTT 


da. 


cosh2y T + sinh2y T cos2a T 
nJ nJ J 


= i 


(2.51) 


This analysis therefore adds a further link to the theory of intensities 
in discrete and continuous spectra [43-44] by including the quasi-bound 
states . 

The above general expressions are readily reduced to Breit-Wigner 

2y j 

form over the sharp resonance region (e nJ >> 1) , by the linear 
expansion 


a.j ~ (n + - )?r + 


3a, 


3E 


(E - E ) 
nJ 


“ < n + I ),r + H5 (E - E nJ> • 


(2.52) 


with the results 


nJ res) (E) = arctan[r nJ /2(E nJ - E) ] 


(2.53) 


T< res) (E/; - 


r .fi 

2Trfi | A(E,J)| 2 = y ~2 ^ 


, r T + (E - E V 
4 nJ nJ 


(2.54) 
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where 


r _ = (2hoj/'iT)e 

nJ 


-2Y. 


nJ 


(2.55) 


Connor [28] and Dickinson [29] would obtain the same result directly 
from (2.41) by applying the outgoing boundary condition 


C" = e^cosctj - ie ^sinot^ = 0 , 


(2.56) 


and the linear expansion (2. 52) to obtain the complex energy of the 
Siegert state. 


E = E nJ - j (2hai/ir)e 


"^nJ 


(2.57) 


-2y. 


nJ 


^ 1 • r nJ is 


Note that in this sharp resonance limit, e 
directly related to the barrier transmission probability given (2.40) 
in the semi-classical form [3]. 


F nj - ra/t vib > (2 - 5 « 

where t = 2 t t/w . 
vxb 

It remains to determine the parameters y, 0 and cj> in (2.39), for 
different barrier models. The simplest theory employs independent 
semi- classical connection formulae at the turning points b and c 
[22-23] with the results 


Y T = /°|k T (r) dr + Jin2 

nJ b 1 nJ 




(2.59) 


0 


(2.60) 
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where the subscript on k nJ (r) implies that the energy is chosen to 
satisfy (2.48). This may be termed the simple semi-classical estimate. 
The first improvements [24, 26, 28] employ a quadratic approximation at 
the barrier maximum 


V T (r) = V wax - jk(r - r ) 2 , (2.61) 

J J z. max 

to obtain 

Y T = cosh _1 [l + e _27T E ] 1/2 , (2.62) 

nJ 

4> nJ = (e) = r( j + ie) - e£n| e| + e , 

(2.63) 

9 nJ = ° * 

where 

e = (E^ - v“ ax )/ha)* 

to* = (7t/m) 1/2 (2.64) 


Methods employed by other workers give the same form for y but 
either a different phase correction [13], or an undetermined phase 
connection [25, 27, 30]. Values of the function c])( 0 which determines 
the level shift in (2.49) and some limiting properties [13, 26, 45] are 
given in table 2.1. A function of this form also plays an important 
part in scattering theory by removing a classical singularity in the 
deflection function [13]. Equations (2 . 62)- (2 . 64) define what we shall 
term the JWKB quadratic approximation. Miller and Good [46] have shown 
how this quadratic approximation may be replaced by a more general 


0 
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association between the physical barrier and the quadratic form (2.61), 
the association being defined by the identity 

e = " ? 'bl k nJ (r) l dr ’ (2 ‘ 65) 

which takes place of (2.64). Our name for this improvement is the 

associated quadratic approximation. First suggested without proof by 

Connor [45], it reduces to (2.64) in the quadratic approximation and 

causes the expression (2.62) for y^. to go over to the simple semi- 

— 2tt£ 

classical form (2.59) when e >> 1 (recall that e < 0 for 
energies below the barrier maximum). Since however [13,27,45] 

<J>(e) ~ e 1 + -jg|o £ 3 for e > 1 » (2.66) 

the level shifts predicted by equations (2.49) and (2.63) (but not 

covered by the simple semi-classical theory) depend inversely, rather 

than exponentially, on e . Hence the leve shift may be orders of 

magnitude larger than the line width in the sharp resonance region. 

Dickinson [29] and Soop [47] have taken the association method one 

stage further in order to obtain a two parameter description. Dickinson 

2 

[29] relies on an association between k T (r) and the inverted Morse 

nJ 

form 

t 2 {q 2 - [2e _ * - e _2s ] } , (2.67) 

2 

where q is the ratio of the available energy to the barrier height 
(measured from the dissociation limit) and t is defined by the 


equation 
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rrt ( 1 - q) = /£|k nJ (r)|dr , (2.68) 

where, as in (2.65), k nJ (r) contains the true potential function. 

The expressions (2. 62)-(2. 65) are then replaced in the present 
notation, by 

y T = sinh ^[e TTE ' - cosh tte /sinh 2irqt] (2.69) 

nJ + 

<j>nj = X(2qt) - <Ke_) (2.70) 

0 nJ = _ <j>(e + ) (2.71) 

where 

e ± = t(l ± q) (2.72) 

X (y ) = argF (iy) - y£n|y| + y + tt/ 4 = x^ 2 y) ~ <Ky) (2.73) 

and (e) is defined by (2.63). This we shall term the associated 
Morse approximation. Note that £_ is identical with -£ defined by 
the quadratic association formula (2.65). Furthermore t is typically 
large unless the barrier is very low ( t lies in the range 4-15 for J 
levels of the ^ ground state considered in table (2.41)), and 
X(y) - -l/12y for y > 2 . Hence it may be verified that (2 . 69)- (2 . 72) 
reduce to (2. 6 2) — ( 2 . 63) with t given by (2.65) except for very low 
barriers (low J values) or for (very sharp) levels close to the base 
of the barrier. 

2 

A final somewhat similar (associated sech ) form may be obtained 

2 

from the results of Soop [47] based on association between k (r) and 
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t 2 {q 2 - sech 2 s} 


(2.74) 


The parameters q and 
Dickinson, but now 


Y 


nJ 


t have the same definitions as those of 


' h _1 / 
smh /- 


cosh[fr(t - 7 


sinhuqt 


— } 


0 nJ = 2x(qt) - <KO + <]>(£_[.) + & 

0=0, (2.75) 

nJ 


where 


- <t‘ - 


i ) 1/z * „ 


6 = tin 


- (t 2 - y ) 1 / 2 in(e'/e') - qtin 
1 - q h "t* - 


e|e^ 


l_ (i - q 2 )t 2 


(2.76) 

§ 2.3 Quasi-bound energy levels 

Any attempt to calculate the width of a quasi-bound state requires 
prior knowledge of its position. Calculated level positions have been 
reported for the ground states the [14, 15, 17, 19, 20, 22, 23] 

HD and D 2 [20], OH [19] 4 HeH + , 3 HeH + , 4 HeD + , 3 HeH + [5], and 
ArK and RbK [31] systems and for a standard model based on the 
Lennard-Jones potential [16, 21, 29]. 

The first calculations [14, 15, 17] employed the points of inflexion 
in the phase shift, r)j(E) > derived by numerical integration of (2.3), 
but this approach has been superceded by use of the time delay function. 
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T d (E,J) [20], after recognition that T^CEjJ) may be derived from 
the same scattering solution, normalised in the form 

r co 

X £J (r) ~ (4m/fik) ^ 2 sin{kr - Jtt/2 + ^(E)} , (2.78) 

by the quadrature [35], 

t d (E,J) = f Q ( j ^ EJ (r) j 2 - I <P E j (°°)/ y^dr + (m/hk 2 )sin(2n J (E) - Jtt) 

(2.79) 

Hence the resonance points, and particularly the level widths, may be 

obtained more conveniently from the maxima in T d (E,J) . The positions 

of peaks in the amplitude ratio, A(E,J) , first demonstrated by 

Buckingham and Fox [14], have also been employed [18,20], Results 

derived from A(E,J) for the ground state of H are found to agree 
-1 

within 0.1 cm with those derived from x d (E,J) [20] for the sharp 
resonances (T < 5 cm ^) but discrepancies, due to variation of the 
ambient phase shift n°(E) in (2.22), of the order of 1-6 cm ^ are 
reported for the broad levels [18, 20] (see table 2.3). Resonance 
positions discussed below refer to the maxima in A(E,J), as having 
the greatest spectroscopic significance. 

A major problem encountered in these calculations is that in seeking 
the scattering solution (such that \p„ T (0) = 0) there is no efficient 
algorithm for converging on a resonance^ a problem which becomes acute 
for very sharp resonances [20] ; the method proposed by Johnson, Balint- 
Kusti and Levine [48] suffers from the same disadvantage. Bain and 
Bardsley [21] have therefore devised a different exact procedure 



23 


specifically for the treatment of sharp resonances. Two solutions of 
(2.3) are introduced, one of which, f (r) , is defined to have zero 
derivative at the origin and to behave as 

f(r) r ^ sin(kr + 6) (2.80) 

at infinity. The second solution, g(r) , is chosen such that 

y ->- CO 

g(r) ~ cos(kr + 6) (2.81) 


and integrated back to the origin. The resonances are then shown to 
occur at energies for which g(0) = 0 , and the level widths for 
rotationless states (J = 0) to be given by an expression analogous to 


(2.79), 

- 2 . 

r _ = Aim 

nJ m _ ^ , 


{[/*( g 2 (r) - f 2 (r))dr - ^ sin(2k7?+ 26) ] -1 } . 

(2.82) 


A similar more complicated expression is given for J ^ 0 . A compari- 
son between the level positions obtained for the Lennard-Jones model and 
those derived from the phase shift is given in table 2.2. 

Other more approximate methods have also been suggested. One 

approach is to convert the problem into a bound state form by imposing 

a boundary condition either at the barrier maximum [19, 20], or at the 

outermost turning point [5, 20]. Studies on the H ^ ground state [20] 

(see table 2.3) favoured the requirement that the wavefunction should 

behave as an Airy function of the second type (increasing exponentially 

into the barrier region) at the outermost turning point; a similar 

suggestion [5] is that ip (r) should have the same logarithmic 

Xj J 
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derivative at this point as the irregular Bessel function Y -^(kr) 

J+ 7 

[41]. This approach still requires numerical solution of the Schrodinger 
equation. The semi-classical theory of J 1 2.2 shows that this may be 
replaced merely by two quadratures (over the potential well a < r < b , 
and the barrier b < r < c ) at each point on a suitable energy grid. 
Dickinson [29] has applied the modified quantisation condition equation 
(2.4), with (0 + cjj) given (2.63) and (2.65), to obtain level positions 
for the Lennard-Jones model in good agreement with the exact values 
(see table 2.2). It is interesting that almost equally good agreement 
may be achieved with a different phase correction by taking the JWKB 
theory to second order [16]. Bayliss [31] has also obtained good agree- 
ment with exact numerical results [49] for the heavier KA.r and RbAr 
systems using this second (now very small) phase correction. 

^2.4 Level widths 

Given the resonance position, computation of the level width 
presents little difficulty. Values have been obtained from the slope of 
Dj(E) [14, 15, 17], from the width and resonance value of x^(E,J) [20], 

from the width [18, 20] and the peak value [5] of the amplitude ratio 
A(E,J), from the special formula of Bain and Bardsley [21] and from 
the semi-classical equations (2.55) and (2 . 59) -(2 . 76) [29, 32]. 
Comparisons are given for the Lennard-Jones system [16, 21, 29] and 
for the ground state of [20] in tables 2.2 and 2.3 respectively. 

The most striking feature is the accuracy of the very rapid JWKB 
estimates. The values obtained by Dickinson [29] for the Lennard-Jones 
system in the associated Morse approximation (equations (2.55), 
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(2.68)— (2. 72)) agree to within 10% with the exact values of Bain 
and Bardsley [21] , and with the probably less accurate values derived 
from the slope of the phase shift by Bernstein et al. [16]. This agree- 
ment is possibly not surprising since these resonances are very sharp. 

Table 2.3 also shows close agreement between the 

exact line widths and the simplest semi- classical estimates (equation 
(2.59)) [20]', Clhd this agreement may be improved for the broader levels 
(r > 20 cm ^) , with no more effort by use of the more refined associated 
parabolic equations (2.62) and (2.65) in place of (2.59). Thus the 
simple semi-classical estimate ma y improved to a value 

r (1) by the algorithm 

r (1) = (2ftto/ir)e~ 2 ^ 

Y = cosh -1 {[l + hw/2nT (0) ] 1/2 } . (2.83) 

Values of Hoo , defined by (2.36), required for the computation of 
r 0) were not reported [20], but estimates may be obtained from the 
tables of accurate level positions given by LeRoy [50]. In fact quite 
rough estimates will suffice, since the level of agreement between 
r (1) and the exact values shown in table 2.4 (IF^^ - F I < 

0.5 cm except for the two broadest levels, (v,J) = (0,38) and (8,21)) 
is unaffected by an increase in hoo by 100 cm 2 . The final columns 
of table 2.4 also contain values of Dickinson's parameters q and t 
defined by (2.68), in order to demonstrate by the magnitude of fc 
( TTt » 1) the essential equivalence between the associated quadratic 
and associated Morse approximations for the levels in question. Hence 
the associated quadratic approximation which requires the three quadratures 
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defined in (2.36), (2.42) and (2.65), may be taken to be the optimum 
semi-classical form for most purposes. 

From an experimental viewpoint the most important theoretical con- 
clusion is the very strong (exponential) energy dependence of the level 
width. Thus, as confirmed by experiment [51], a typical rotational 
progression for a hydride system will terminate in absorption in the 
form sharp, broad, very diffuse, followed by a complete loss of structure 
which may obscure the existence of one or possibly two quasi-bound 
(E < V” ax ) upper state levels. This must result in a discrepancy 

between the curve of limiting predissociation and the locus of the 

THciX 

barrier maximum, Vj , which is estimated by LeRoy and Bernstein [20] 

to range from zero at J = 0 to as much as 1000 cm ^ at J = 38 for 
the ground state of H 2 . Hence although extrapolation of the limiting 
curve should yield the correct dissociation limit [6,10], deductions 
about the long range behaviour of the potential [52-54] must be treated 
with caution. Any attempt to extrapolate from observed data to the 
true locus of V™ ax requires an understanding of both the measurably 
sharp and the very broad levels discussed above. 

§ 2. 5 Summary and conclusions 

The properties of non-physical Siegert states with complex energy 
levels, exponentially decreasing time evolution and purely outgoing 
characteristics have been used in 2.1 to unify the resonance behaviour 
of the phase shift the collision delay time and the spectroscopically 
observed internal to external amplitude ratio. The existence of such 
states also explains the Breit-Wigner parameterisations of n T (E), 
r d (E,J) and A(E,J) (equations (2 . 21)- (2 . 24) ) . 


c 
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A general semi-classical theory leads first inf 2. 2 to a quantisation 

condition (equation (2.48)) and to expressions for the levelshift, (2.49), 

and for the line width (2.55) in terms of three parameters y T , Q r 

nJ nJ 

and <f> n j . Formulae for these parameters are then given for four different 
treatments of the barrier region, of which the associated quadratic model, 
equations (2.62), (2.63) and (2.65), appears to be the most generally 
applicable. 

Different techniques for the location of quasi-bound states are 
discussed inf 2. 3. The most significant practical conclusion is that the 
(associated quadratic) semi-classical quantisation condition appears to 
be remarkably accurate. 

Recent semi-classical formulae for the level width are also shown 
to achieve high accuracy even for very broad levels. Thus there would 
appear to be little necessity for exact numerical computations in this 
field. The existence of reliable analytical formulae may also be expected 
to lead to more accurate location of barrier maxima, in cases where the 
highest quasi-bound levels are indiscernably diffuse. 

III. PREDISSOCIATION BY INTERNAL EXCITATION 

Just as predissociation by rotation may be seen as another 
manifestation of a shape or orbiting resonance, so is predissociation 
by internal (electronic) excitation (Herzberg's type I [6]) associated 
with the Fano [69] or Feshbach [70] resonance phenomenon in scattering 
theory. The situation is complicated by the necessity to consider at 
least two electronic states. Physical mechanisms and selection rules 
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governing the necessary mixing of electronic states have been well 
characterised [8, 9] for many years. Our present interest is in the 
origin of the detailed dependence of the line width or lifetime 
pattern as a function of vibrational and rotational state. Recent 
work [56-68] has confirmed an early suggestion by Rice [55] that 
fluctuations in the line width, due to interference between the nuclear 
wavefunctions for the interacting state, provide a sensitive measure 
of forms of the relevant potential energy curves. 

The subsequent theory is set in context by a short resume of the 
Fano [69] theory of interaction between discrete levels and the 
continuum; alternative approaches have been suggested by Feshbach [70] 
and more recently by Van Santen [71]. We also include for the sake of 
completeness a brief account of possible electronic origins of the 
interaction term. 

3. 1 The Fano lineshape 

Fano [69] first considers the coupling between a single discrete 

level and a continuum of states ip . Second order interaction 

n E 

between <j> n and other discrete states , are therefore neglected 

at this stage. The normalisations are chosen such that 

(<J> j H | <j> ) = E , 

n 1 1 n n 

(ip E , |H| Tp E „) = E'6(E" - E') 

(V |H|4> n ) = (<J> n |H|i]> E ,)* = V E , n . -(3.1) 


This implies according to (2.5) that 
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t 


E' 




2m 


E' 


TffiV 


1/2 


sin[k'r + 6 (E f ) ] 


(3.2) 


where iii_ t depends on the other (electronic and angular) coordinates 
E 

of the system. Our purpose is to find an (energy-dependent) 
eigenfunction 


V E = a(E)<j) n + / b(E J E , )^ E ,dE' , (3.3) 

subject to the same asymptotic normalisation. The functions a(E) and 
b(E,E') will then determine the line shape. 

It follows from (3.1) and (3.3) that at the energy eigenvalue E , 


E^a(E) + j V nE ,b(E s E')dE' = Ea (E) (3.4) 

v w , a(E) + E’b (E,E' ) = Eb(E,E') . (3.5) 

iii n 

b(E,E') is now eliminated in favour of an admixture ratio Z(E) 
which becomes the central quantity in the theory, by a rearrangement 
of (3.5) [36]; 


b(E,E’) 


E ~ E , + Z(E)5(E - E’) "| V E , n a(E) . 


(3.6) 


At the same time the presence of the resulting pole due to the term 
(E - E') ■*" makes it necessary to specify how the integrations in (3.3) 
and (3.4) are to be carried out; otherwise Z(E) is not uniquely 
defined. The most convenient choice is the Cauchy principal value [72], 
denoted by P below. 

Hence on substitution for b(E,E') in (3.4) and removal of the 


common factor a(E) 
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where 


Z(E) 



F(E) 




(3.7) 


(3.8) 


It remains to determine a(E) by applying the asymptotic normalisation 
condition (3.2) to in (3.3). Equations (3.2), (3.3) and (3.6), 

£j 

together with the identity 


p r°° sin(k'r + 6) 
(E - E') 


dE' 


P / 


2k’ ) si 

k + k'J 


sin[(k - k')r + kr + 6] , 

(k - k 1 ) 


P / 


f 2k' ) 


k + k' 


,, , sin(k' - k)r „ , 

cos (kr + 6) (k k 1 ) ^ 


-tt cos(kr + 6) , 


(3.9) 


imply that 


XT -»• co 


% 


2m 


-i 1/2 


„ 2. 
Tffi k 


a(E)V„ { — 7T cos(kr + 6) + Z(E)sin(kr +6)} 
bn 


where 




2m 


•7 1/2 


TTh 2 k 


a(E)V [tt 2 + Z 2 (E)] 1/2 sin(kr + 6 + 6 (r) ) , 
bn 

(3.10) 


tan 6 (r) = — Ctt/Z (E) ) , 


(3.11) 


XT co 

because (p 0 . Hence for the correct asymptotic behaviour 

n 


a(E) = vj [TT 2 + Z 2 (E)] 1/2 


(3.12) 
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This function, which represents the amplitude of the discrete component 

<p in the energy normalised solution (3.3), is the direct analogue of 
1 n 

the amplitude ratio A(E) of the previous section 

With Z (E) given by (3.7), equations (3.11) and (3.12) show that 
6 and a(E) have the characteristic Breit-Wigner resonance behaviour s 
given by (2.21) and (2.22), 




arctan[F /2(E' - E) ] , 
n n 


(3.13) 


a(E) 


2tt 


n 


(E - 


2 12 
E') + t r 

n 4 n 


(3.14) 


where 

E^ = E n + F (E) (3.15) 

F n ■ 2T I V J 2 • 

Furthermore by analogy with (2.51) 

r E’+r /2 

n n | a(E) j dE = 1 . (3.17) 

/ 2 

n n 

Notice that the level width depends, via V^ E , only on the 

resonance continuum state ip , but that the levelshift function F(E) 

c* 

appears in (3.9) as an integral over the full continuum. Both T and 

n 

F(E) are strictly energy dependent but they may be evaluated with 

sufficient accuracy at the resonance point E if the line width is 

n 

small compared with the level spacing. 
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Finally Fano [69] obtains a compact expression for the spectro- 
scopic transition moment by collecting the above results in the form 

\ = ^n Sin C0S * (3.18) 

nE 

where represents the original discrete state <j> modified by an 

admixture of continuum states, 

V nF' V 

$ = $ + p/ dE' . (3.19) 

n n E - E 

Hence the transition moment between and another state 6 becomes 

E r o 

^e^IV = sin 6(r> - ^ e |t U 0 > cos 6(r) ■ 

nE 

(3.20) 


The spectroscopic absorption line shape therefore differs according 

to the relative values of (0 I T | d> ) and (U;_ I T I d> ) . It takes the 

n o E o 

Lorentzian form 


l<f E l T !Vl 2 


r 

n 

2ir [ (E - E') 2 + T 2 /4] 
n n 


«J t M 2 


(3.21) 


when transitions between <p and the continuum if; are forbidden, but 

o E 

appears as a Lorentzian emmission line, against the continuum absorp- 
tion, if ($ I T | <j> ) = 0 ; 
n 1 1 o 


(¥ e |t| 4> o ) 



r 2 /4 

n 


[(E - E’) 2 + r 2 /4] 
n n 




(3.22) 

Interference between the two terms in (3.20), results in a distorted 
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line shape (see figure 3.1) characterised by Fano [69] in terms of the 
ratio 

q = ( $ n | T | <J> o ) / ttV* e ( | T | 4> o ) . (3.23) 

Fano a -t' 

Herzberg and Jungen [73] have observed this characteristic/Beutler [74] 

line shape in the pre-ionised Rydberg spectrum of • 

Fano [69] also extends the argument to cover several discrete states 

4> n interacting with one continuum ^ , and several continua 

interacting with one discrete state. The former naturally always 

applies in diatomic systems, but corrections to the one level analysis 

above become significant only when the line width given by (3.16) becomes 

comparable with the level spacing. The latter situation may arise in 

certain diatomic molecules, but it is of more general significance in 

the polyatomic case because each (vibrational-rotation) channel has an 

associated continuum of states [75]. 

The major correction to the theory in the many level, one continuum 

case is that the perturbed levels appear as the eigenvalues, , of 

a level shift matrix [69] with elements Ed + F (E) , where 

n nm nm 


F 

nm 


(E) 


P / 


V nE’Vm dE ' 
E - E' 


(3.24) 


This leads to the introduction of transformed discrete state combinations 


*v - i %A < 3 - 25 > 

where the c form the eigenvectors of the above matrix. The energy 
normalised solution finally appears in the form 
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= cos 


5< r >{£ 


Ev 


V E - E 




V 


(3.26) 


where 


r(r) _ y ^ V Ev 
tan o = l 


V E - E. 


V 


. . V vE'V dE ' 

% ‘ + p/ — e--" r - 


(3.27) 

(3.28) 


The conditions for the validity of (3. 15)- (3. 18) are therefore 

first that |F (E) I << E - E , so that 
1 nm 1 n m 


E - E + F (E) , (3.29) 

V n nn 

and secondly than jv^j <<c - in which case the sums in (3.27) 

and (3.28) may be approximated by a single term at each resonance, and 

the term [cos v / (E - E )] reduces to (sin 6^ ; /irV ) as in 

.EiV v nb 

(3.18). 


The results for the one discrete, several continua case may be 
summarised in the form [69] 


$ sin 6 (r) 
™nE n 


v (i) 

r nE 

■ V 
i nE 


,(i) 


cos 


(r) 


(3.30) 


where 


= <j> + £ P 


n 


V nE'»E' )dE ' 
E - E' 


nE 


= Z V 


(i) | 2 

nE 1 


(3.31) 


and 


is given by (3.13). Both the level shift F(E) 


and the 
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line width therefore appear, via V ^ in (3.25), as a simple 

sum of the previous single continuum formulae taken over the continua in 

question. 

§ 3.2 The interaction term 

The theory in the previous section is quite general. When specialised 
to a diatomic system the interaction term may be written 

V nE' = o X 2n^ H 21^ X 1E^ dr 5 (3.32) 

where 

H 21 (r) = <P 2 |H(q,r)|P 1 > . (3.33) 

Here n denotes a particular vibrational, total angular momentum level, 
|P 1 > and |P 2 > are used to designate appropriate electronic and 
rotational angular momentum states^and the integral in (3.33) is taken 
over the electronic coordinates q and the angular coordinates of the 
nuclear position variable r . X 2n ( r ) and X 1E ( r ) are the corresponding 
vibrational and continuum wavefunctions , normalised to unity and to a 
delta function of energy respectively. 

Predissociation may be induced by both internal (intramolecular) 
and external interactions. The important internal cases are covered 
by the selection rules and expressions for H-^Cr) given by Kronig 
[8]. These require in every case 

AJ = 0 , +*-*+ , g«-*u , (3.34) 


but interactions governed by 
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AA = 0,±1 in cases (a) and (b) 

A(2 = 0,±1 in case (c) (3.35) 

are allowed. There is also a weak selection rule 

AS = 0 in cases (a) and (b) 

which may be broken by spin-orbit coupling. An important observational 
distinction is made between homogeneous and heterogeneous inter- 
actions [76] for which AA = 0 (or AQ = 0 in case (c)) and AA = ±1 
(or Afl = ±1) respectively, in that the latter^which arise from rotational- 
electronic coupling^ yield an expression for H^ 2 ( r ) proportional to 

1/9 

[J(J + 1) - A(A ±1)] .A short table of expressions for 
evaluated in the approximation of pure precession has been given by 
Czarny, Felenbok and Lefebvre-Brion [61]. 

Observed predissociations have also been attributed to the external 
interactions induced by magnetic fields [9, 77] and molecular collisions 
[80-82]. A condensed form of Van Vleck's theory [9] of magnetically 
induced predissociation has recently been given by Chapman and Bunker 
[79]. Qualitative aspects of the theory of pressure induced predissociation 
have been considered [81, 83] but a comprehensive theory remains to be 
developed. 

^3.3 The Franck-Condon approximation 

The existence of an intersection between the repulsive and bound 
state curves, V^r) and V 2 (r) respectively, may be used to justify a 
significant simplification in equation (3.32), because the effective 
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integration region is known [84] to be localised at the crossing 
point, R; an analytical estimate obtained in ^3. 5 suggests an effective 
range between ±0.05 £ and ±0.20 X in typical systems. Under these 
conditions it may be permissible to replace H^ 2 (r) by its value, H ° 2 > 
at the crossing point, so that 

r n = 2TT{V nE* 2 = 2^|H° 2 | 2 |/“x 2n (r) X 1 E (r)dr | 2 . (3.36) 

The validity of this Franck-Condon approximation to the interaction 
element is supported by exact numerical calculations [56,61], and by 
analytical arguments in §3. 5 below. The use of a similar Franck-Condon 
approximation in the computation of the level shift by means of 
equation (3.8) is however open to greater risk because the significant 
integration range in the evaluation of V^, varies as E' spans the 
continuum. 

A second limitation on the validity of (3.36) is that the line 
width obtained should be small compared with the level spacing (see 3.1). 
A significant case not covered by (3.36) therefore arises when H^ 2 (r) 
is sufficiently large that the predissociation occurs from the upper 
adiabatic potential curve, V + (r) , where 

V±(r) = \ [V x (r) + V 2 (r)] ± j {[V^r) - V 2 (r)] 2 +4H 2 2 (r)} l/2 , 

(3.37) 

rather from the diabatic curve V 2 (r) . The observed predissociation of 
ICl [85] and IBr [86, 87] are assumed to be of this type. In these 
circumstances equation (3.16) remains valid provided that x 2n ( r ) and 
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Xie ^ r ) i- n (3-32) are replaced by the bound, y | ^ (r) , and continuum, 
X „(r) , wavef unctions determined by V (r) and V (r) respectively 
and the interaction term is written as [60, 88] 


H + _(r) = -H_ + (r) 


d 2 e de d_ 

dr dr dr 9 


(3.38) 


where 0(r) is the adiabatic mixing parameter 


0(r) = arctan 


2H 12 (r) 
V x (r) - V 2 (r) 


However the presence of a peak in (d0/dr) , and the disappearance of a 
crossing point precludes the use of a Franck-Condon approximation in the 
evaluation of V , given by (3.32). 

ilk 

Intermediate cases between these near diabatic and near adiabatic 
limits may be covered by diagonalisation of the Fano level shift matrix 
defined by (3.24), in order to determine the appropriate admixture (3.25) 
of bound states. A simpler alternative may be to seek a direct 
solution of the appropriate coupled differential equations, either by 
numerical [60] or by analytical [66] techniques. It is interesting to 
find that an exact model computation of this type, covering the full 
range of possible interaction strengths [60], shows very few resonance 
positions which depart significantly from these obtained by either a 
strict diabatic (distortion) or adiabatic approximation. This study 
gave however only a graphical indication of the level widths. 

Finally attention may be drawn to the possibility of predissociation 
in the absence of a potential curve crossing (type a° , b° or c° in 
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Mulliken's classification [7]). Such cases may fall into either the 
"non-diabatic" (predissociation from ^^(r)) or "non-adiabatic" 
(predissociation from V (r)) category, and no decision can be made 
without examination of possible changes in the electronic state. Even 
in the non-diabatic case however the Franck-Condon approximation (3.36) 
would require careful justification in view of the absence of a crossing 
point . 

3.4 Exact numerical results in the Franck-Condon approximation 

Within the validity of the Franck-Condon approximation (3.36), a 
computation of the line width rests on knowledge of the exact numerical 
wavefunctions X lE ^ r ) and X 2n ( r ) • Henceforth n will be taken to 
include a pair (v,J) of vibrational rotational quantum numbers. 

There is no problem such as that encountered for shape resonances (see 
2.3) in locating the level position because x 2n ( r ) may be taken as a 
vibrational-rotational eigenfunction of the experimental RKR bound 
state curve ^ 2 ( r ) . For any given repulsive curve V^(r) , the 
continuum state X]_ E ^ r ) as a l so r eadily determined; the only problem 
concerns normalisation. The use of a true continuum state normalised 
as in (3.2) by comparison with a spherical Bessel function in the 
asymptotic region [89] is strictly preferable on grounds both of compu- 
tional speed and flexibility, but a pseudo-continuum approximation is 
frequently employed [58,59], according to which 

X lN (r ' ) ’ 


X iE (r ) 


3N 


3E. 


N 


(3.40) 
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where X-^ N ( r ) is the N-th normalised bound wavef unction for the 

potential V^(r) modified by an infinite wall at some large distance. 

Given such wavef unctions , it is possible to vary the assumed repulsive 

curve until the computed line width pattern agrees with that obtained 

by experiment. Computations of this type have been performed for 

3 _ 

predissociations from the states of 02(B [56, 57], N 2 (b ' ^ ) 

[62], H 2 and D 2 (D'tt ) [58], 0D(A 2 Z + ) [61] and I 2 (B 3 tt( 0+) [58, 63]. 
A similar calculation for the levelshift has been reported for the 
B(0*) of Se 2 [64]. 

The first general conclusion is that the results are indeed sensi- 
tive to the form of the repulsive curve. This is drmatically illustrated 
by the calculations of Reiss and Ben Aryah [57] and Murrell and Taylor 
[56]. Secondly the line width pattern is found, as first predicted by 
Rice [55], to fluctuate with v at a frequency which is high or low 
according to whether the intersection occurs on the attractive or repul- 
sive branch of V 2 (r) (types c + or c in Mullikeh's notation [7]). 

This behaviour may be attributed to a variation in the relative phases 
of the wavefunctions and ^2n^ r ^ at t * ie crossing point. Hence 

the oscillation frequency depends on the rate of divergence of V^(r) 
and V 2 (r) above the crossing point. Typical patterns obtained by 
Murrell and Taylor [56) are shown in figure 3.2. The irregular nature 
of the fluctuations in case c + means that a single repulsive curve 
may give rise to two or more regions of strong predissociation separated 
by intermediate v levels for which the predissociation is weak. 


3 _ 

Theory [56, 63] suggests that the observed predissociations of C> 2 ( £ ) 
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3 + 

[90-91] and I 9 [ tt( 0 )] [63, 78] may be understood in this way. A 

^ u 

second consequence is that even when the predissociation becomes so 
strong that all levels appear indiscernibly broad one may expect to 
find a few scattered sharp lines; the visible spectrum of IBr [87] 
is a clear example of this. The observation of such "vestigial remains" 
[55] would itself provide important information about the repulsive 
curve . 

Similar fluctuations in the theoretical line width as a function 

of the rotational quantum number (in this case N) have been reported 
2 + 

for the A( £ ) of OD [61]. The variation with N(N + 1) is now 

however sufficiently slow that a more regular oscillatory pattern may 

be recognised (see figure 3.3). Again, within the accuracy of the 

measured line widths, the experimental pattern is found to be consis- 

• 4 " 

tent with a single repulsive curve of type c , attributable on the 

•k 4 

evidence of the magnitude of line widths to the state OTTO ( tt) . A 

similar calculations for the heterogeneous b predissociation from 

the D ( ’ TT ) to the B'('I + ) of H. [59], shows over the small J range 
u n 2 

considered (0 < J < 10) no Franck-Condon contribution to the line 
width variation with J ; the entire variation obtained is may be 
attributed to the J dependence of H °2 in (3.36). 

The above conclusions are derived in the main for predissociation 

+ + + 
of types c . For an extension to cases a and b“ for which the 

asymptotes lie at and above the curve crossing point respectively [7], 

one may note that the predissociation pattern depends by the above 

arguments only on the local forms of the potential curves in the 
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crossing region. Hence one must expect an abrupt onset of the fluctuating 
line width pattern in place of the smooth increase shown in figure 3.2. 
Furthermore any consistent theory must also account for any perturbations 
below the predissociation limit. No calculated predissociations of 
types a 1 , b 1 and c 1 , for which the crossing point coincides with the 
minimum in the bound state curve [7], have been reported but there ap- 
pears no reason to expect a qualitatively different line width pattern. 

For the reasons given above, no conclusions can yet be reached for the 
cases designated a° , b° and c° . 

Turning to the level shift, Atabek and Lefebvre [64] have recently 

applied equation (3.8) to the interpretation of reported shifts in the 

3 _ 

B( state of Se2 [93] . The experimental results are shown to be 
consistent with an intersection of type c + , but a discrepancies are 
found between these Fano level shifts, and the Born Oppenheimer (or 
adiabatic) shifts attributable to the difference between ^(r) and 
the lower adiabatic curve V_(r) defined by equation (3.37). The 
significance of these discrepancies remains to be investigated. Other 
features of interest are first that the calculated level shifts below 
the crossing point (|F(E )| - 8,10 cm ^ for v = 18,19) exceed the 
corresponding widths (T^g < 0*00 cm \ - 0.02 cm ^) by several 

orders of magnitude. Secondly a projection to levels above the crossing 
point indicates that the level shifts oscillate about zero as v 
increases. Quantitative conclusions drawn from a calculation of this 
type must be regarded as tentative in view of the lack of any objective 
criterion for the definition of an experimental level shift. Nevertheless, 
it is tempting to suppose that simultaneous predictions about the 
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behaviour of the line width and the level shift above the crossing 
point might lead to experimental detection of the fragmentary remains 
of the spectrum which could be used to refine the details of the model, 
even when most lines are indiscernibly diffuse. 

Analytical approximations which shed some light on the above 
numerical conclusions are developed in the following section. 

§ 3.5 Analytical approximations 

Analytical approximations to the line widths and level shifts are 
conveniently developed in the notation indicated in figure 3.4. Here 
V d (r) denote the centrifugally corrected potential curves 

V^(r) = V ± (r) + (J + ) 2 -K 2 /2mr 2 , (3.41) 

1 2 

in which the Langer substitution [40] of (J + — ) for J(J + 1) , 

(made here for consistency with the discussion in 2.3) may be regarded 

2 2 

as a common small semi- classical correction, -ft /8mr , to the two 
potential curves. This substitution has no effect on the position, R , 
of the crossing point, with energy 

E XJ = V 1 (R) + < J + | ) 2 -ft 2 /2mR 2 = V 2 (R) + (J + | ) 2 fi 2 /2mR 2 , 

(3.42) 


or on the potential derivative difference 


where 


AF = F - F 

1J 2J * 


F iJ " -< 3 Y dr >r=R ’ 


(3.43) 


(3.44) 



44 


which play an important part in the theory. AF is in fact independent 
of J , and (in the notation of figure 3.4) necessarily positive. 
Finally it is convenient to define two phase integral combinations 


<J ‘ k U <r)dr + 1 2J k 2J (r)dr 

XJ K 


<f> T = / R k (r)dr - / R . , . , 

vJ a 2J 2J a U k U (r)dr 


(3.45) 


for the interpretation of + and - predissociation respectively; 
here the functions k.^.(r) and k 2 ^(r) are defined by (2.2) and (2.36) 
with E = in each case. By this definition <j>^j may be associated 

with the limiting form of the upper adiabatic potential defined by (3.37), 


V d (r) « V d (r) r < R 

* V d (r) r > R 


(3.46) 


as H 12 (r) -*■ 0 . 

Franck-Condon approximations . The first analytical line width expressions 
are obtained in the Franck-Condon approximation, applicable for small 
H °2 in (3.36). The simplest is that derived by the method of Landau 
and Lifshitz [84] from the semi-classical forms (2.30) and (2.32); thus 


2mKoj 


X 2vJ (r) 


vJ 


Tffi 2 k 2J (r) 


1/2 


sin[/ 

l r 


2J 


k 2J (r)dr + tt/4] 


W r) ' 


2m 


Tm"k 2J (r) 


1/2 


sin [/ r k (r)dr + tt/4] (3.47) 


‘ a U 1J 
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+ + *4* 

for a predissociation of type a , b or c . The argument involves 
a quadratic expansion about the crossing point for the phase of the low 
frequency component of the product ^ 2 v j ( r ) in (3.36), leading 

to the following form for the line width 


vJ 


^yjKzl 2 

im 2 v 2 

vJ 


cos(0x 2 + <j)^j + tt/ 2) dx 


= sin 2 (*+ 3 + tt/4) , 


(3.48) 


where 


vJ 


-ftv AF 
vJ 


(3.49) 


and 


B = AF/2Rv 


vJ 


1 * 
T mV T 

2 vJ 


E T - E t 
vJ xJ 


(3.50) 


A similar expression 


r vJ = C sin2( Kj + 


(3.51) 


is obtained in cases a , b or c 

The relation between the Fresnel integral in (3.48) and the error 
function [41] indicates that 95% of the integral arises from the range 

i i 1/2 

|x| < 2(2/8) . Hence equations (3.48) and (3.51) depend on the 

validity of the semi-classical approximations (3.47), and on the assump- 
tions of linearity in V^(r) and constancy in over the range 
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r = R ± 2[ttv vJ /AF] 1/2 

~ R i 5[(E vJ - E xJ )/mAF 2 ] 1/4 , (3.52) 

where in the second line of (3.52), the units are taken to be r(£) , 

E(cm "*■) , m(atomic mass units), F(cm ^ £ 2 ) . The significance of 

this transition zone, the width of which increases with increasing 

energy, was first recognised by Bates [94]. Parameters applicable to 

3 - 3 + 2 + 

the predissociations of 0^( E ) , ^(B 17 and OH(A Z ) are 

listed in table 3.1. Note that since AF may become quite small for 
inner (a , b or c ) curve crossings, the Franck-Condon approximation 
may require special justification in these cases, particularly for 
systems of low reduced mass. 

A second analytical approximation is available to cover energies 
near and below the crossing point, for which overlap between the 
transition zone and the classical turning point regions invalidates the 
semi-classical approximation (3.47). One is necessarily now interested 
only in predissociation of types c + and c . In the latter case for 
example, the following Airy function [41] approximations are applicable: 


*2vJ (r> “ 


2mRto 


vJ 


a 2J* 


1/2 


Ai [ -a 2j(r 


a 2J )] 


X lEJ (r) 


2m 


“li* 


1/2 


Ai[-aij( r - a u )) , 


(3.53) 


where 


a u " (“ij /* 2 ) 173 


(3.54) 
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and the normalisations are chosen to be consistent with those in 
(3.47) (see Landau and Lifshitz [95] ). Equations (3.53) define exact 

wavefunctions for the potential curves 


V. _(r) = E — F (r — R) 

xJ xJ lJ 


(3.55) 


with F > 0 . Hence the turning points at energy E are given by 
iJ vJ 


a. T = (E . - E _)/F. t . 
iJ vJ xJ iJ 


(3.56) 


With the necessary overlap integral derived in the appendix, the final 
line width expression may be shown to be 


r _ = Trr° T (v T /v*)Ai 2 [-(E T - E J/E J , 

vJ vJ vJ J vJ xJ J 


(3.57) 


where 


E J = \ mV J 2 = ^ 2 Fij F 2j / 2iiiAf2)1/3 • (3.58) 


This formula was first given by Rice [55], and rediscovered independently 

by Degenkolb , Steinfeld, Wasserman and Klemperer [78] and the present 

author [65]. It may be shown to go over to (3.51) when ( F v j “ F x j) >:> 

•k 

Ej because [41] 


Ai(-Z) Z ~ > 1 tt 1/2 Z 1/4 sin( | Z 3/2 + tt/4) 


(3.59) 


and within the validity of (3.55) 



^ 4 ^ U R [E T 
* a 2J vj 


- / R [E T - E . 
a^j vJ xJ 


E xJ + F 2J< r - R)]1/2dr 
+ F 1J (r - R) ] 1/2 dr 


2 

3 


* 3/2 

I(E vJ - E xJ )/E j' • 


(3.60) 
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Equations (3.57), (3.58) and (3.60) (with in place of ) 

may also be shown to be valid in case c ; the argument requires merely 

a reversal of sign in the argument of X 2 v j( r ) and th e introduction 

of l F 2jl ;i ' n t * le definition of a^j . 

As a measure of the relative validities of (3.48) or (3.51) and 

(3.57) we may note that equation (3.59) reproduces the Airy function 

within ±0.025 at the first maximum, Ai (-1.012) = 0.1536 , and within 

±0.005 for Z < -3.0 . Hence as a general rule the semi-classical 

forms (3.48) or (3.51), which require a linear approximation to the 

potential only over the transition zone (see (3.52)), may be taken to 

* 

be superior for - E^j > 1.5 Ej . 

Before examining the physical significance of these results, we 
may note that an extension of the above argument to include an 
exponential interaction term 


H]_ 2 (r) = ‘ H i 2 exp[-a(r - R) ] , 


(3.61) 


may be shown to yield [96]. 


r = Trr 

VJ vJ 


/ . ' 
V -r 


'fi 2 a 3 (F, T + F„J 


- 

E T - E T 

vJ 

exp 

v 1J 2J' 

Ai 2 


vJ xJ 

V* 

vJj 

3m AF 3 


e j 


- E r 


12 


(3.62) 


where 


E^ = (fi 2 a 2 |F u F 2J |2niAF 2 ) . 


(3.63) 
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The asymptotic form of 




exp 


(3.62), 


fl2a3<F U + F 2J> 


- 

r \ 

E , - E 

. 2 
sin 

2 

vJ xJ 

3m AF 3 

3 



3/2 


2 -« 

a nv . 


AF 



(3.64) 


differs from that derived from (3.57) by the presence of the exponential 

2 

factor and the additional phase term, a -fiv T /AF . The former is 

V vJ 

2 3 3 -3 

probably unimportant in most cases since ["K a (F.j + F2j)/3mAF ] = 10 
for F^j = 1600 cm ^ , F^j. = 8000 cm , a = 1 A and m = 10 a.m.u 

for example. The form of the phase correction again underlines the 
significance of equation (3.52); it will be small provided the range, 
a ■*" , of H^^) is large compared to the width of the transition zone. 

We turn now to consider the physical factors affecting the line 
width variation with v at given J . Figure 3.5 gives a comparison 
between the predictions of equations (3.48) and (3.57) and the computed 
values of Murrell and Taylor [56] . The properties of the Airy 
function [41] indicate that the line width envelope should reach a 

JU 

maximum at E^j + 1.019 , with a point of inflexion at the crossing 

point and an exponential decrease 


i> r vjl ex P 



[E - E 

3/2" 

4 

xJ vJ 


3 

E* 



J J 



(3.65) 


for E - E^j < -Ej . The much slower decline at higher energies is 
governed according to equation (3.49) by the ratio (w /v ) ; note 

VJ v J 
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that since the local energy spacing, > must vanish at the 

dissociation limit, even when the predissociation appears very strong^a 
discrete spectrum must be regained as the energy increases. 

The frequency of oscillations within this envelope depends on the 
energy variation of the appropriate phase integral or <f> v j • Seen 

in this light the oscillations may be attributed to semi-classical inter- 
ference [97] between the two possible escape trajectories associated 
with transitions at the crossing point during outward or inward motion 
respectively, the latter being followed by reflection at the classical 
turning point on the repulsive curve. Hence the line width pattern is 
in effect an interferogram for the nett path difference b^j - a^ or 
a lJ ~ a 2J cases + or - respectively. This idea is fully 
exploited in£3 . 6 . 

Seen in more concrete terms, it is clear from figure 3.4 and from 

i 

the definitions (3.45) that the energy dependence of (J)“ must depend 
on the rate of divergence of the potential curves. In case c + the 
necessary derivative may in fact be associated by the usual semi- 
classical arguments (see (2.36)) with a hypothetical energy spacing 
hu/ + ^ in the potential curve V^(r) defined by (3.46); 


SE _ 3<j> _ 

vJ vJ 


3v 3E 




(3.66) 


This means that close to the dissociation limit, where both hco and 

vJ 

ho/ + ^ must depend only on the long range form of V"?(r) [98], 


= 7T . 


(3.67) 
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The factors affecting the line width variation with J may also 
be understood in the light of equations (3.48), (3.51) and (3.57). The 
overriding consideration is that the energy difference in (3.57) varies 
as 

E t -E t = E- E + [hcB - -R 2 /2mR 2 ]J(J + 1) . (3.68) 

vJ xJ v xo v 


Hence for vibrational levels near the crossing point the line width must 
increase or decrease with increasing J in cases c + or c 
respectively. Similarly the phase terms and <j> may be seen to 

be increasing and decreasing functions respectively, with (9<J> + /9.J) 

VJ 

approximately given by [67] 


3<j). 


+ 

vJ 


3J 


= (2J + 1)tt 


B - B (+) (E ) 

V V 




(+) 


J 


(3.69) 


where hu/ + ^ and B + (E^) are respectively the hypothetical vibrational 
spacing and rotational constant at the measured energy E^ for the 
potential curve V + (r) . 

Finally the effects of isotopic substitution may be considered. 

It is convenient to make comparisons at a common absolute energy and at 

equivalent mass reduced J values (having a common value of 

J(J + l)/m) . With these provisos, the envelope function in (3.48) 

1/2 

and (3.51) carry however a factor m . Hence when regarded as a 
function of energy the onset of strong predissociation will appear 
sharper and the subsequent oscillation frequency will be increased at 
higher mass values. 
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General approximations . A quite different type of semi-classical theory, 
analogous to that described in 2.2 has also been developed [66]. This 
relies on the use of semi-classical connection formulae (equation (18) 
of reference [66]) for the Stuckelberg-Landau-Zener [99-101] curve 
crossing model. The results are therefore again limited by the approxi- 
mations of constant velocity, constant interaction and linear potential 
curve expansions over the transition zone defined by (3.52). Being 
non-perturbative, however, the method carries no restriction on the 
magnitude of • Furthermore it yields analytical expressions both 

for the line width and for the level shift. 

The results may be expressed in terms of one or other of the 
Landau-Zener transition probabilities [99-101] 

P _ = 1 - exp[-2TTH° ' 2 /Kv AF] 

vJ 12 vJ 

P V J - 1 - P vJ • exp[-2ra° 2 2 /*v vJ AF] : (3.70) 

P^j gives the probability of a transition from V^r) to V^(r) on a 
single passage through the crossing zone, while P^ refers to a 
transition from one adiabatic curve V^(r) or V^(r) to the other. 
Predissociations of measurable line width are predicted when either 


P _ « 1 or P _ « 1 . 

vJ vJ 


(3.71) 


The former correspond to what we have termed the non-adiabatic pre- 
dissociations discussed in the first part of this section, and formulae 
for the line widths are identical with those given by (3.48) and (3.51). 
The associated level shifts may be written 



53 


AE 


l r O . , ,+ , IT > ,,+ . IT > 

= 2 r vJ sln ^ vJ + 4 > cos ^ vJ + 4 > 


vJ 


(3.72) 


in cases 


+ 


+ + 

b or c 


and 


AE VJ = “ \ r vJ sin( Kj + I } COS (< ^vJ + J ) (3 - 73) 

in cases a , b or c 

The second inequality in (3.71) refers to predissociation from the 
upper adiabatic potential curve V d (r) defined by (3.37), such as that 
observed in the visible spectra of IC1 [85] and IBr [86-87]. It is 
convenient to coin the notations a , b or c" to distinguish predis- 
sociations of this type from the more common non-adiabatic phenomena. 
The expressions for the line width and leve now become 


r r o 2 ,(2) 

r T = T cos (j> T 
vJ vJ vJ 


AE ° \ r° Sin (j) (3) cos <j> (2) 
vJ 2 vJ T vJ vJ 


(3.74) 

(3.75) 


where 


few 


vJ 


= 2 


vJ 


exp[-2iTH° 2 2 /hv vJ AF] , 


(3.76) 


and 


P (2) 

V J 


r b 2J , , . , (2m) <L , b 2J r _ „J. . 1 1/2, 

= L, k 2J (r)dr = — L [E vJ- V 2 (r)] dr • (3 ‘ 77) 


2J 


2J 


Note that "R w v j now refers to the level spacing in V (r) . The ratio 

of the line width to the local energy spacing, r° /-ftw , is therefore 

VJ VJ 

predicted to oscillate within an increasing envelope as indicated 


schematically in figure 3.6. The disappearance of hw 


vJ 


at the 
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dissociation limit however again implies that the lines must 
eventually become sharp. 

These equations raise several points of interest. First the 

inequalities (3.68) give quantitative tests for the validity of a 

non-diabatic or a non-adiabatic description. Secondly equations (3.74) 

and (3.77) predict a sharp line in the spectrum whenever a bound level 

of V^(r) coincides with a bound level of V^(r) because, by the Bohr 

( 2 ) 1 

quantisation condition (2.34), <J> = (n + -r )tt under these conditions. 

VJ ^ 

This behaviour is strikingly demonstrated by the visible spectrum of 
IBr [86, 87]. Finally the intimate connection between the line width 
and the level shift which implies zero level shift for the sharp 
resonances under all conditions may have some general significance. 
Although this behaviour is not confirmed by the numerical results of 
Atabek and Lefebore [64] possibly due to strong curvature in the repul- 
sive curve V^(r) , it means within the validity of the present model 
that the measured energy variation with J contains a systematic error 
in the rotational constant. Thus at the sharp resonance points 


3E r 
vJ 

3J(J + 1) 


hcB ± | r° 
v 2 vJ 


3ij>; 


vJ 


3J(J + 1) J 


(3.78) 


for non-diabatic predissociations, while 


9E 


vJ 


3j (J + 1) 


= hcB r° 
v 2 vJ 


3<J>. 


( 2 ) 

vJ 


9J (J + 1) J 


(3.79) 


in the non-adiabatic case. The behaviour of [ 3 ^ v j- / 3 J ( J + 1)] has been 


discussed above (see 3.69)); the application of similar arguments shows 


that 


e 
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B - B (2 ^(E ) 
v v 

( 2 ) 

Iko 


This term will be negligible in cases a , b , and c because the 
minima in V+(r) and V^(r) must coincide, while [3<{/ 2 V8J(J + 1)] 
will clearly be negative in cases a , b , and c . Hence equations 
(3.69) and (3. 78)-(3. 80) carry a general prediction that anomalies of 
this type must lead to an increase in the apparent B value. The 
magnitude of this increase clearly depends on F°^ or F°^ . It must 
be negligible for predissociations which are weak to the extent that 
all rotational lines are observed, but it may become significant in the 
interpretation of fragmentary systems, such as that shown by IBr [86, 
87], which contains only a few measurably sharp lines. 


§ 3.6 Direct inversion of the line width pattern 
A direct method for determination of the repulsive curve V^(r) 
and the interaction strength H° 2 has been developed [65, 67, 68] from 
the analytical theory of the previous section. Closely akin to the 
RKR method [102, 103] of bound state spectroscopy, it is based on 
recognition that the phase integrals $ ^ in equations (3.45), (3.48) 
and (3.51) contain the same information about the upper branches of 
V^(r) and V^(r) in figure (3.4), as that supplied by the Bohr quanti- 
sation condition (<j> = (n + ■— )tt) for a normal bound state curve. The 
only difference is that the energy dependence of must be extracted 

from the line width pattern. Once this function is known however the 
separation between the turning on V^(r) and V 2 (r) at energy U is 
given by [65, 68] 
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b 2 (U) - a 1 (U) 


+ 


in case c or [68] 


a x (U) - a 2 (U) 


- 1 

' 2 \ 1//2 r.u 

S4> + 

vo 

dE 

TT ! 

mj J E 

xo 

l SE J 

(U - E) 1/2 

■K 

f l] 1/2 r-u 

r s<r ] 

vo 

dE 

TT 

1-J J E 

XO 

l 3E J 

(U - E) 1/2 

turning point 

a 2 (U) or b 2 (U) 


(3.81) 


(3.82) 


known from an experimental RKR curve V^(r) , equations (3.81) or (3.82) 

are sufficient to determine V^(r) . Similarly the J dependence of 
i 

4> v j may be used [67] to obtain an analogue of the second RKR equation, 
but this is unlikely to be of practical value except for confirmatory 
purposes. 

+ 

The determination of <|)“ rests primarily on location of the 
crossing point E xo . As a first step it is convenient to rewrite 
equations (3.48), (3.51) and (3.57) in the form 


Y = (r /hu> ) 

VO VO VO 


1,2 - tt 1/2 KE*- 1/4 Ai 

O 


E - E 

VO XO 


E* 

o 


for E - E 

VO XO 


(3.83) 


K(E - E ) 1/4 sin(c]> + y + n tt) 

vo xo vo 4 v 


for E » E 

vo xo 


(3.84) 
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where 


oj 

vo 


3G 

v 


3v 


(3.85) 


Y^ o and E are therefore experimental quantities; the unknowns are 

* ± * 

K, E , E , tp and the integers n . Two equations for K, E 
o xo vo ° v o 

and E^ o may be taken from the position and magnitude of the first 
Airy function maximum estimated from three neighbouring line widths (or 
from the greatest measured with in the absence of adequate data [68]); 
the equations are 

max = 7r 1/2 Ai (_ 1-01879 ) K E * _l/4 = 0.94943K E* _1/4 (3.86) 

o o 


E max _ E = 1.01879E* . 

xo o 


(3.87) 


For a third equation it may be assumed that the maximum product 

Y vo ( E vo “ R x0 ) 1/4 in the higher energy region corresponds with a 

sinusoidal maximum (a crude estimate of E will suffice for this 

xo 

purpose) . Hence 


K = [(E - E ) 1/4 y ] 

vo xo vo max 


(3.88) 


As an alternative, depending on the quality of the data one might employ 

four measured line widths in the Airy function region to determine the 

first maximum and a point of inflexion in the line width envelope, the 

latter being taken to determine the crossing point. Preliminary estimates 
* 

°f K > E^j and E^. obtained in this way may be improved by comparison 
between computed and measured line width patterns in the Airy region. 

The phases in this region are now given by equation (3.60); 
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2 * * 3/2 

4 [<e - e )/e y u 

3 vo xo o 


(3.89) 


It remains to eliminate ambiguities of n^TT in the phases obtained 
by inversion of (3.84). A plot of 

+ ^2 + 

Ad) = d> - / k.(r)dr in case c 

T vo vo 2 

K 

R 

= / k„(r)dr - $ in case c (3.90 

2 vo 

3/2 

against “ E q ) is recommen ^ e ^ f° r this purpose. Here k^(r) 

is derived from the known RKR curve ^(r) by means of equations (2.2) 
and (2.6). This graph is suggested by a linear expansion for the 
potential V^(r) (see equation (3.60)); the present model requires that 
it should pass through and behave linearly at the origin. Not only will 
this graph remove the necessary ambiguities of n^TT , it may also be 
used to smooth out experimental uncertainties in the line width, to 
refine the location of the crossing point in cases offering insufficient 
data in the Airy region, and to combine data for different isotopic 


substituents because the phases in (3.90) all vary at a given energy as 

1/2 


Finally one may note an internal check on the consistency of the 

* 

method, because the value of E q used to fit the data in the Airy 

function region, may be compared with that deduced from (3.58) in 

terms of the slopes F. and F„ derived from the calculated curves. 

lo Zo 

A valuable general property of the method is its stability. Refin em ent 

of a preliminary analysis of the magnetic flourescence quenching data 
3 -{- 

for I 2 [B f( 0 u )] [63, 78], by means of the A<J> plot (figure 3.7) is found to 

* _1 

change the crossing region parameters (E , E and K) from (306 cm . 

xo o 5 



59 


-1 -1 -1 
468 cm and 3.217) to (475 cm , 340 cm and 3.162) , but to alter 

the turning points on V^(r) by -.002 X for v = 6,11 , ±.001 £ fo 
v = 14,16,18,21,24,28 and 32 and -.004 £ for v = 38 . There is in- 
sufficient experimental data near the predicted crossing point for a 

JU 

meaningful consistency check on the value of E . A comparison between 
the refined potential curve V^(r) an ^ t ^ ie ^ orm deduced by Chapman and 
Bunker [63] from a numerical Franck-Condon calculation is given in 
figure 3.8. 

This direct inversion procedure has been successfully applied to 

3 - 

the computed predissociation data for 0 ^( £ ) [65] and to experimental 

2 + 3 + 

data for 0D(A £ ) [67] and (B tt ( 0^) ) [68]. The resulting repulsive 

potential curves are closely comparable to those deduced by numerical 
Franck-Condon calculations (see figure 3.8). The advantages of the 
direct approach are speed and flexibility, since the only prior assump- 
tion is of linearity in the potential curves over the transition zone. 
This suggests its general future use as a preliminary to the potentially 
more exact numerical Franck-Condon procedure. 

§3.7 Summary and conclusions 

It has been amply confirmed by practical tests that a pattern of 
measured predissociation line widths (or lifetimes) contains sufficient 
information to determine the repulsive curve V^(r), over the predis- 
sociation region, and also the interaction potential • The line- 

shape envelope as illustrated in figures 3.2, 3.5 and 3.6 suffices to 
determine the predissociation type. The frequency of oscillations 
within this envelope depends, in the normal non-diabatic case, on the 
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rate of divergence of the two potential curves above the crossing point. 

Similar fluctuations in the case of predissociation from the upper 

adiabatic curve V + (r) depend on the shape of V£(r) • Subsidiary 

fluctuations with J are also predicted, at a frequency roughly 

governed according to equations (3.69) and (3.80) by the displacement 

between the crossing point and the bound state potential minimum. An 

analytical characterisation of these fluctuations has led to development 

of the direct inversion procedure described in 3.6. The interpretation 

+ 

of a well developed predissociation pattern in cases c is therefore 
now a routine matter. 

Preliminary attention has also been directed towards the concurrent 
levqlshift. Specifically, at energies below the crossing point one may 
expect to observe shifts which are several magnitudes larger than the 
line width. However it is not yet clear whether these should be inter- 
preted as Fano shifts in the sense of equation (3.8), or simply as 
adiabatic shifts due to an avoided curve crossing. The sign of such 
shifts is negative for crossings of the + type, and may be expected 
by the adiabatic interpretation to be positive in the opposite case. 

At energies above the crossing point the levelshift is predicted to 
oscillate about a mean value zero in such a way that within the validity 
of a Landau-Zener model a small positive additional term is predicted in 
the rotation constant. 

Turning to the future one may expect to see the direct connection 
between perturbations and predissociations of types a and b more 
fully exploited, particularly when the whole spectrum is substantially 
discrete. A more challenging situation for the theoretician arises 
when the predissociation is strong in the spectroscopic sense that 
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all discrete structure appears lost, but weak to the extent that the 
line width is small compared with the vibrational spacing. Here one 
might hope that analysis of low-energy level-shifts or possible pertur- 
bations might lead, via the appropriate repulsive curve, to the pre- 
diction and detection of fragmentary remains of the spectrum at the 
oscillation minima in figures 3.5 and 3.6. 

Another possibility is that spectroscopic predictions of strong 
predissociation into states which correlate with ground state atoms 
might stimulate the observation of such resonances by scattering 
techniques, as a complement to the achievements of Bernstein [1966], 
Stwalley et al. [3], and Schulte et al [4] in the somewhat simpler 
shape resonance field. A possible example might be the IBr system 
[86, 87], but there may be others which are preferable on experimental 
grounds . 
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APPENDIX: THE AIRY FUNCTION OVERLAP INTEGRAL 


The integral 

I = /“ Ai[y(r - c)] Ai[y'(r - c’)]dr , (A.l) 


with y > y' may be evaluated by use of the integral representation 
[41] 

1 i , 3 

Ai(z) = — / cos( -r u + uz)du = -x— / exp[ — h iuz]du . 

iro 3 2 tt r 3 

(A. 2) 


Thus 


I 


1 

( 2tt ) 2 


d d d ex p{ 3 ^ + + i(Yu + y'v)r 

i(uyc + vy' c' )}drdudv 


exp{ <u 3 + v 3 ) - i(uyc + vy'c')} 6(yu + y'v)dudv 


2^7 d ex P^ f t 1 “ (Y'/Y) 3 ]v 3 - i(c’ - c)y'v}dv 


(Y 3 - Y' 3 ) 1/3 Ai 


YY'(c - c T ) 


, 3 ,3 N l/3 

(Y - Y ) 


(A. 3) 


Here the identity [36] 


CO 

/ exp ikr dr 
—00 * 


2tt 6 (k) 


(A. 4) 


has been used to obtain the second line of (A. 3). 



68 


CAPTIONS FOR TABLES 

2.1 Values of the level shift function <j>(e) in the associated 
quadratic approximation (equation (2.63)). 

<j>(-e) = -<K £ ) • 

<Ke) - e _1 /24 + 7e" 3 /2880 for £ > 1 . 

2.2 Comparative resonance positions and widths for the Lennard-Jones 
model, with the potential parameters employed by Bernstein et al. 

[6]. £ q denotes the well depth. Energies and widths are expressed 
as a fraction of the well depth. 

2.3 Comparative resonance positions and widths for the ground 

state [10]. Values under BC(Airy) refer to the Airy function 
boundary condition at the outermost turning point. 

2.4 Simple, J and improved, , semi- classical level widths 

for the ground state, q and t in the final columns are 

the parameters for the associated Morse approximation; note that 
Trqt » 1 . 

3.1 Parameters for some observed predissociations. 
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CAPTIONS FOR FIGURES 

2.1 The condition for possible predissociation by rotation. 

2.2 Notation for the semi-classical theory, a, b and c denote the 
classical turning points. Capital letters denote the amplitudes of 
travelling waves, with primes and double primes used to indicate 
outgoing and incoming motion respectively. 

3.1 The Beutler-Fano line-shape for q = 3 (see equation (3.23)). The 
corresponding form for q = -3 is obtained by reflection in the 
origin. 

3.2 Computed linewidth variations with v [57] for cases c + (upper 
diagram) and c (lower diagram). Note the contrast between the 
irregular fluctuations in case c + , with the smooth variation in 
case c . Compare also figure 3.4. 

3.3 Computed linewidth variations with v and J [61] for homogeneous 

2 + 

predissociation from the A( £ ) of OH . 

3.4 Notation for the semi-classical theory (a) in case c , and (b) in 
cas e c 

3.5 A comparison between the computed linewidth pattern [57] for pre- 

3 - 

dissociation of 0^( E^) and the analytical expressions (3.48) 
(dashed line) and (3.57) (solid line). Note that a satisfactory 
Airy function representation (solid line) as far as the third 
minimum should be regarded as atypical. Curvature in the potential 


\ 
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energy functions may normally be expected to cause deviations between 
the semi-classical (dashed line) and the Airy function behaviour at 
a point between the first maximum and the first minimum (see text). 

3.6 Schematic line width variation with v for a c + case of non- 
adiabatic predissociation from the potential curve V (r) . Vertical 
lines indicate the line widths. 

3.7 Refinement of the magnetic flourescence quenching data for the 

3 + 

B[ tt(C> u ) ] state of I . The point v = 6 is derived a postiori 

* 

from equation (3.60) with E q given in terms of the slope of the 
inverted curve by equation (3.58). 

3.8 Comparison between the repulsive curve V^(r) deduced for the 

3 + 

B Tr (° u ) state of by a numerical Franck-Condon calculation [63] 

and the points obtained by the direct inversion method. 
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TABLE 2.1 



e 

<j>(e) 

£ 



0 

0 

0.6 

0.080 

0.1 

0.137 

0.7 

0.068 

0.2 

0.150 

0.8 

0.058 

0.3 

0.135 

0.9 

0.051 

0.4 

0.115 

1.0 

0.045 


0.5 


0.096 



i 


E nJ 



Phase Shift 

BB^ 

JWKB^ 

118 

.31005 

.31004 

.31003 

197 

.31005 

.31008 

.31007 

25 

.34855 

.34862 

.34846 

77 

.34935 

.34935 

.34938 

128 

.35005 

.35007 

.35007 

123 

.38967 

.38961 

.38962 

87 

.39006 

.39009 

.39008 

137 

>.4012 

.40233 

.4025 


Bernstein et al. [6] . 
Bain and Bards ley [11]. 
^ Dickinson [20], 




Phase Shift BB^ JWKB (c) 





1.2 

X 

10 17 

1.0 

X 

10 17 




9.6 

X 

10~ 28 

9.0 

X 

io“ 28 

6.2 

X 

IO " 4 

6.0 

X 

io" 4 

6.7 

X 

10~ 4 




6.9 

X 

O 

1 

oo 

8.0 

X 

10 " 8 




2.0 

X 

10 -11 

2.0 

X 

io" 11 

1.4 

X 

io" 4 

8.4 

X 

io -5 

1.0 

X 

1 

o 

rH 

4.1 

X 

io -4 

4.2 

X 

10~ 4 

4.1 

X 

10 4 




3.4 

X 

io -3 

3.2 

X 

io" 3 





TABLE 2.3 


V 

J 


E nJ <cirl > 


T T (cm 
nJ 

" ; S 

T , (max) 
d 

A (max) 

BC(Airy) 

A(E S J) 

•JivcCB 

0 

38 

7510.0 

7514.0 

7508. 7 

98.1 

8 V . 0 

0 

37 

6513.3 

6513.5 

6513.3 

5.98 

5 . 60 

1 

35 

5549.8 

5550.0 

5549.7 

.14.3 

14.9 

2 

33 

4688.4 

4689.0 

4688.2 

20.8 

2 '/. ■ 5 

3 

31 

3925.0 

3925.4 

3924.9 

25.1 

26. 1 

4 

29 

3254.7 

3255.4 

3254.8 

25.4 

28.4 

5 

27 

2673.0 

2673.8 

2673.4 

25.8 

29.3 

6 

25 

2175.0 

2176.0 

2175.7 

27.4 

31.4 

7 

23 

1755.3 

1756.7 

1756.4 

31.7 

36 9 

8 

21 

1407.0 

1409.4 

1409.4 

42.1 

4 G . 3 

9 

19 

1121.6 

1127.2 

a 

66.2 

a 

9 

18 

725.9 

725.9 

726.0 

0.53 

0.55 

10 

16 

586.0 

586.0 

586.1 

2.93 

3.22 

11 

14 

480.1 

481.0 

481.7 

18.5 

22.3 

11 

13 

199.4 

199.4 

199.4 

0.005 

0.0053 

12 

12 

385.0 

398.6 

a 

11.6 

tr.-. 

12 

11 

215.5 

215.5 

215.6 

2.62 

3.09 


-.7 
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TABLE 2.3 (cont'd.) 


V 

J 


E nJ (an ’ 1) 



X 

• (max) 
d 

A(max) 

BC(Airy) 

A(E s J ) 

JWKB 

13 

9 

195.0 

205.6 

a 

89.6 

a 

13 

8 

89.9 

90.0 

i— 1 

o 

CTk 

1.89 

2.38 

14 

6 

81.9 

121.0 

a 

79.0 

a 

14 

5 

45.7 

49.2 

a 

26.4 

a 

14 

4 

3.76 

3.76 

3.76 

0.007 

0.0085 

£ 

The boundary 

condition places this level above 

the barrier 

maximum . 
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TABLE 2.4 







. 



V 

J 

■hoj 3 

r (o)b 

r (l)c 

r b 

exact 

d 

q 

t 


0 

38 

1050 

87.0 

69.8 

98.1 

.9983 

62.5 

l 

35 

1000 

14.9 

14.2 

14.3 

.9 75.’ 

15.2 

2 

33 

950 

22.5 

21.0 

20.8 

.9775 

13.5 

3 

31 

900 

26.7 

24.5 

25.1 

.9 778 

12.1 

4 

29 

875 

28.4 

25.8 

25.4 

,9773 

11.1 

5 

27 

800 

29.3 

26.3 

25.8 

.9765 

9.97 

6 

25 

750 

31.4 

27.8 

27.4 

.9767 

9.13 

7 

23 

675 

36.9 

31.7 

31. 7 

.9 799 

8.44 

8 

21 

600 

48.3 

39.0 

42.1 

.9875 

8. 66 

11 

14 

400 

22.3 

19.1 

18.5 

.9608 

4.27 


Estimates of "Ro) derived by graphical extrapolation from the 
tables of LeRoy [40] have been rounded to ±25 cm The 
level of agreement between and F is unaffected by a 

6X3. Cl 

change of 100 cm ^ in "ftco. 
b LeRoy and Bernstein [10]. 

C Equations (2.83). 


d 


From energy levels in table 2.3 and V given by LeRoy [50]. 
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TABLE 3.1 


State 

Type 

E - E 

X 

AF 

R 

Range 

References 



(cm 1 ) 

(cm ^ 

/Ov 

(A) 

(A> 


°2 

+ 

c 

5000 

50,000 

1.89 

±0.10 

[56, 65] 


c 

3000 

9,000 

2.87 

±0.15 

[63, 68] 


OH A 2 Z + 


+ 

c 


4000 


24,000 


1.65 ± 0 . 20 [61, 67] 
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